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Abstract. The neutrino burst from a galactic supernova can help determine the neutrino mass hierarchy and 613, and provide 
crucial information about supernova astrophysics. Here we review our current understanding of the neutrino burst, flavor 
conversions of these neutrinos, and model independent signatures of various neutrino mixing scenarios. 
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Neutrino flavor conversions inside a SN are sensitive 
to extremely small On values and the nature of neutrino 
mass hierarchy 1 1]. The observation of the neutrino burst 
from a galactic SN will therefore provide information 
complementary to that from a long baseline experiment. 
It will shed light on many of the outstanding questions in 
neutrino oscillation physics and astrophysics. 

In this article, we shall follow the production, prop- 
agation and detection of SN neutrinos, which are akin 
to the source, the long baseline and the far detector of a 
neutrino factory setup. 

OPERATION OF THE SN v-FACTORY 
Core collapse and SN explosion 

Neutrinos and antineutrinos of all species are pro- 
duced inside the SN through pair production processes. 
In addition, v e is also produced by electron capture on 
protons: pe~ — > nv e . Before the collapse, neutrinos of 
all species are trapped inside their respective "neutri- 
nospheres" around p ~ 10 10 g/cc. 

When the iron core reaches a mass close to its Chan- 
drasekhar limit, it becomes gravitationally unstable and 
collapses. A hydrodynamic shock is formed when the 
matter reaches nuclear density and becomes incompress- 
ible. When the shock wave passes through the v e neu- 
trinosphere, a short v e "neutronization" burst is emit- 
ted, which lasts for ~10 ms. The object below the shock 
wave, the "protoneutron star," then cools down with the 
emission of neutrinos of all species. This emission takes 
place over a time period of t ~ 10 s |2[]. 

The eventual explosion of the star involves the stalling 
of the original shock wave, its revival by the trapped 
neutrinos, and a "delayed" explosion where large scale 



convections play an important role Jj, |4|] . However, for 
the ^10 sec neutrino burst that we focus on here, the 
actual explosion mechanism, or even whether the star 
successfully explodes or not, is mostly immaterial. 



The source: primary neutrino fluxes 

A SN core acts essentially like a neutrino blackbody 
source with flavor-dependent fluxes. Since the fluxes are 
almost identical for V^V^Y^ and v T , all these species 
may be represented by v T . The "primary fluxes" F® may 
be parametrized by the total number fluxes OofVa), av- 
erage energies (Eo(va)), and the "pinching parameters" 
that characterize their spectral shapes |5J]. 

The values of the parameters are highly model de- 
pendent, as can be seen from Table [T] The two leading 
models - the Livermore simulation [6] and the more re- 
cent Garching calculation [7] - agree on (Eo(v e )) « 12 
MeV and (Eo(v e )) « 15 MeV, and have consistent val- 
ues for the pinching parameters, but they differ widely on 
(Eq(v x )) and the ratios of total fluxes. In particular, the 
equipartition of energy assumed in the Livermore model 
is not a feature of the Garching model. 



TABLE 1. Differences in flux predictions 
from SN models 



Model 


(Eo(v x )} 


<t>o(Vv) 


4>a(v e ) 
4>o(y,) 


Garching 
Livermore 


18 
24 


0.8 
2.0 


0.8 
1.6 



In the light of the model dependence, it is important to 
make sure that the inferences drawn from the observed 
neutrino spectra do not depend strongly on the exact 
model parameters. 
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FLAVOR TRANSFORMATIONS ALONG 
THE LONG BASELINE 

Neutrinos that are produced approximately as mass 
eigenstates at very high ambient matter density in the 
core propagate outwards from the neutrinosphere. They 
have to travel through the core, mantle and envelope 
of the star, through the interstellar space, and possibly 
even through some part of the Earth before arriving at 
the detector. Inside the SN, collective and MSW mat- 
ter effects take place. In the interstellar space, neutrino 
mass eigenstates travel independently, whereas oscilla- 
tions, enhanced by MSW effects, occur inside the Earth. 



MSW resonances inside the SN 

The traditional picture of flavor conversions in a SN 
is based on the assumption that the effect of neutrino- 
neutrino interactions is small. In this case, flavor con- 
versions occur most efficiently in the MSW resonance 
regions. SN neutrinos must pass through two resonance 
layers: the H-resonance layer at Ph ~ 10 3 g/cc char- 
acterized by (A/WatjQ, On), an d me L-resonance layer at 
Pl ~ 10 g/cc characterized by (Am 2 , 612). This hierar- 
chy of the resonance densities, along with their relatively 
small widths, allows the transitions in the two resonance 
layers to be considered independently J^] . 

The outcoming incoherent mixture of vacuum mass 
eigenstates is observed at a detector as a combination of 
primary fluxes of the three neutrino flavors: 

Fv e = pF° e + (l-p)Fl, (1) 
F% = PFI + {1-P)FI, (2) 

where p and ~p are the survival probabilities of v e and v e 
respectively. 



TABLE 2. Survival probabilities for neutrinos, p, 
and antineutrinos, ~p, in various mixing scenarios 
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The neutrino survival probabilities are governed by the 
adiabaticities of the resonances traversed, which are di- 
rectly connected to the neutrino mixing scheme. In par- 
ticular, whereas the L-resonance is always adiabatic and 
appears only in the neutrino channel, the adiabaticity of 
the H-resonance depends on the value of 613, and the res- 
onance shows up in the neutrino (antineutrino) channel 
for a normal (inverted) mass hierarchy. Table|2]shows the 
survival probabilities in various mixing scenarios. For 
intermediate values of 613, i.e. 10~ 5 < sin 2 9u < 10~ 3 , 



the survival probabilities depend on energy as well as the 
details of the SN density profile |Q]]. 

Scenarios A, B and C are the ones that can in principle 
be distinguished through the observation of a SN neu- 
trino burst. Note that the sensitivity to 613 is comparable 
to the expected reach of a neutrino factory. 



Collective effects at large v densities 

The neutrino and antineutrino densities near the neu- 
trinosphere are extremely high (10 30 ~ 35 per cm 3 ), which 
make the v — V interactions in this region significant 
Indeed, the Hamiltonian is now given by 

fl(p,r) =H vac (j>) + V(r)+H vv (p,r) , (3) 

where H vac is the vacuum Hamiltonian, V (r) is the MSW 
potential due to electrons and the V — V interaction poten- 
tial H vv is OH! [LI 

ffvv(P,r) = V2G F [ 7^3 K-pq (« v p-«vP) • (4) 

J (271) 

Here n v (q,r,f) and p(q,r,f) are the number density and 
density matrix of neutrinos with momentum q, whereas 
K pq = (1 — cos Opq). The quantities with a "bar" represent 
antineutrinos. Note that the evolutions for v and v are 
nonlinear, and are coupled to each other. 

The distinctive features in the flavor evolution of such 
a relativistic gas have been identified in IhI [l5l [l6l \vh . 
The evolutions of the density matrices of neutrinos and 
antineutrinos may be represented in terms of the preces- 
sions of the corresponding Bloch vectors, termed as "po- 
larization vectors" P. The traditional MSW oscillations 
correspond to v and v of each energy independently pre- 
cessing about B, the Bloch vector corresponding to the 
Hamiltonian H vac (p) + V (r). 

Different collective effects play important roles in dif- 
ferent regions of the star [18]. When the neutrino den- 
sity is extremely high, P's of all energies remain tightly 
bound together and precess with a common frequency, 
giving rise to synchronized oscillations 11911 . At lower 
densities, the P's remain bound together to a large extent, 
but have a tendency to relax to the state that has the low- 
est energy, causing bipolar oscillations 1I20I, I2H1 . Collec- 
tive interactions also predict "spectral split," a complete 
swapping of the energy spectra of two neutrino flavors 
above or below a critical energy, as the neutrinos tran- 
sit from a region where collective effects dominate to a 
region where the neutrino density is low B22[ 12311 . 

The analytic treatment of collective effects till now 
has been mostly in the two-flavor limit, assuming a 
steady-state, spherical, half-isotropic, finite source. The 
dependence of the flavor evolution on the direction of 



propagation of the neutrino may give rise to direction- 
dependent evolution J§, QUl, or to decoherence effects 
Halm, but for a realistic asymmetry between neutrino 
and antineutrino fluxes, such effects are likely to be small 
lfl8l l25ll and a so-called "single-angle" approximation 
can be used. Recently, a formalism for analyzing the 
three-flavor effects has also been developed 12611 . 

For SN density profiles where the collective effects are 
over before the MSW effects begin, the collective effects 
are equivalent to changing the primary spectra available 
for further propagation, so that results in Table [2] stay 
valid. For "shallow" density profiles A27I1 where the col- 
lective and MSW regions may overlap, the situation is 
more complex and has to be analyzed separately. 



Oscillations inside the Earth matter 

If the neutrinos travel through the mantle, and possibly 
core, of the Earth before reaching the detector, the neutri- 
nos undergo oscillations inside the Earth and the survival 
probabilities change j28i I29I l30ll . This change however 
occurs only in those scenarios in Table|2]where the value 
of the survival probability is nonzero. This provides the 
means for discriminating between various mixing sce- 
narios. 

When antineutrinos, for example, travel through the 
mantle and the core of the earth, the sharp density jumps 
give rise to the survival probability of the form Bill 

7 

J5 D « cos 2 12 + £A,- sin 2 (<k/2) (5) 

i=l 

in the two-layer model of the Earth, where the coeffi- 
cients A, are functions of the mixing angle 612 in vacuum, 
mantle and core. The phases 0, depend on the distance L 
traveled through the Earth matter and the values of mass 
squared differences Am 2 between Vi and V2 in the man- 
tle and the core. If the neutrinos traverse only through the 
mantle, only one oscillating term is present in (O, with 
<j> = 2Am 2 mantl eL/E JH]. 

DISTINGUISHING BETWEEN 
NEUTRINO MIXING SCENARIOS 



This process, dominant at any water Cherenkov or scin- 
tillation detector, will be instrumental in determining the 
V e spectrum. In order to measure the v e spectrum cleanly, 
however, one needs a large liquid Ar detector, with the 
relevant process v e + 40 Ar — > 40 K* + e~ . 

The uncertainties in the primary neutrino spectra make 
the task of determining the survival probabilities p and 
~p almost impossible, and alternative model independent 
signatures of various neutrino mixing scenarios need to 
be looked for. The Earth effects and shock wave im- 
prints seem particularly promising for this purpose. Since 
the characteristics of the neutronization burst are robust 
across models, the structure of the neutronization peak 
can also identify scenario A, where the peak is highly 
suppressed [34]. 

Earth effects manifest themselves in two ways. Firstly, 
the total number of events and the spectral shape 
changes. Secondly, Earth effect oscillations are intro- 
duced, which may be identified even at a single detector. 

Comparing signals at multiple detectors 

If neutrinos travel different distances through the Earth 
before reaching two detectors, the difference in the sig- 
nals at the detectors could show evidence for Earth ef- 
fects. In order to obtain a statistically significant differ- 
ence in the neutrino spectra, the detectors need to be of 
the size of SK or larger. 

The IceCube detector, though designed to detect in- 
dividual neutrinos with E > 150 GeV, is able to de- 
tect a SN neutrino burst during which the number of 
Cherenkov photons detected by the optical modules 
would increase much beyond the background fluctua- 
tions i35l l36ll . Though this does not measure energies 
of individual neutrinos, the total luminosity can be deter- 
mined at the per cent level. Comparison of the luminosi- 
ties as functions of time at the IceCube and SK (or its 
larger version) can identify the earth effects, since they 
are typically time dependent l3~7ll . 

The relative locations of SK and IceCube imply that 
for the SN in a large portion of the sky, neutrinos 
pass through the Earth for only one of the detectors. 
This makes the SK-IceCube comparison an interesting 
prospect. 



The only SN observed in neutrinos till now, SN1987A, 
yielded only ^20 events. Though it confirmed our un- 
derstanding of the SN cooling mechanism, the number 
of events was too small to say anything concrete about 
neutrino mixing (see [33] and references therein). On the 
other hand, If a SN explodes in our galaxy at 10 kpc from 
the Earth, we expect ~ 10 4 events at Super-Kamiokande 
(SK) through the inverse beta decay process V e p — > ne + . 



Identifying Earth effects at a single detector 

The oscillating terms sin 2 (0,/2) in <(3j can manifest 
themselves as peaks in the Fourier power spectrum of 
the "inverse energy" spectrum of v e 113211 : 
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where y = (25 MeW)/E. The positions of these peaks are 
independent of the primary neutrino spectra, being deter- 
mined by the solar oscillation parameters, the Earth mat- 
ter density, and the position of the SN in the sky. There- 
fore, Earth effects can be identified merely by identifying 
the presence of these oscillation frequencies in the ob- 
served spectrum. If the neutrinos pass only through the 
mantle, there is only one Fourier peak. When the core is 
also involved, three out of the seven possible peaks may 
have significant power ||3lll . leading to an easier identifi- 
cation of the Earth matter effects. 

Finite energy resolutions of detectors tend to smear 
out the modulations in the energy spectrum, and sup- 
press high-A; peaks. The comparison between a simulated 
megaton water Cherenkov detector and a 32 kt scintil- 
lation detector lf3~lTl shows that the better resolution of 
the scintillator detector almost compensates for the much 
larger water Cherenkov detector size. 

The observation of Earth effects in v e (v e ), either 
through the luminosity comparison or through Fourier 
peaks, eliminates scenario A (B) independently of SN 
models. 



NEUTRINOS FOR SN ASTROPHYSICS 
Pointing to the SN in advance 

Since neutrinos are expected to arrive hours before the 
optical signal from the SN, the neutrino burst serves as an 
early warning IB8I1 to the astronomy community. Being 
able to determine the position of the SN in the sky is 
also crucial for determining the Earth crossing path for 
the neutrinos in the absence of the SN observation in the 
electromagnetic spectrum. 

A SN may be located through the directionality of 
the ve~ — > Ve~ elastic scattering events in a water 
Cherenkov detector such as SK 11391, 14011 . The direction- 
ality of this reaction is primarily limited by the angular 
resolution of the detector, the kinematical deviation of 
the final-state electron direction from the initial neutrino, 
and the nearly isotropic v e p — ► ne + background which is 
30-40 times larger than the signal. Adding to the water 
a small amount of gadolinium, an efficient neutron ab- 
sorber, would allow one to detect the neutrons and thus 
to tag the inverse beta reactions 14111 . Efficient neutron 
tagging can improve the pointing accuracy at SK from 
- 8° to - 3° for a SN at 10 kpc & 



Tracking the shock wave in neutrinos 

The passage of the shock wave through the H- 
resonance (p ~ 10 3 g/cc) a few seconds after the core 



bounce may break adiabaticity, thereby modifying the 
spectral features of the observable neutrino flux ll43ll44l 
H^BSl. One expects a "dip" in (E e ) as well as the num- 
ber of events, and a simultaneous peak in (E 2 ) / (E e ) 2 . If a 
reverse shock is also present, the above features become 
a a double-dip and a double-peak respectively 14711 . 

Since the density of the H-resonance layer depends on 
energy, the positions of the dips in the number of events 
at different neutrino energies would allow one to trace 
the shock propagation while it is in the mantle around 
densities of p ~ 10 3 g/cc l47ll . 

The shock wave effects can be diluted by stochas- 
tic density fluctuations as well as turbulence. For ex- 
ample, for 0n > 10~ 4 the shock wave imprints may be 
part ly erased with 5-correlated stochastic fluctuations 
114811 - If the turbulent convection generated behind the 
shock wave is sufficiently large, flavor depolarization 
takes place at the H resonance, so that the sharp shock 
wave effects are replaced by gradual depolarization ef- 
fects H49TI . A recent hydrodynamic simulation 115011 sug- 
gests that some of the shock wave effects survive in 
spite of the smearing factors above. Sterile neutrinos may 
leave their imprints in the shock wave GlL which can 
also survive in the presence of turbulence B52I1 . 

The nonmonotonic density profile of the shock wave 
may cause the neutrinos to pass through multiple H 
resonances. This gives rise to oscillations in the survival 
probabilities of neutrinos, where the positions of maxima 
and minima are independent of primary fluxes 15311 . The 
oscillations are however smeared out due to the finite 
energy resolutions of detectors, and the signal may be 
detectable only in extremely optimistic cases. 

The observation of any of the shock wave imprints in 
the v e (v e ) spectra would imply that the neutrino mixing 
scenario is A (B). 



SYNERGY BETWEEN SN NEUTRINOS 
AND A NEUTRINO FACTORY 

A galactic SN burst is a rare phenomenon, expected to 
occur only once in a few decades. However, its obser- 
vation is expected to reap a rich scientific harvest. It is 
therefore imperative that we are ready with suitable long 
term detectors that will observe the relevant signals. In 
the meanwhile, better theoretical understanding of neu- 
trino transport inside the SN, combined with more ac- 
curate measurements of the neutrino mixing parameters, 
will equip us for making the most of the cosmic catastro- 
phe. 

Earth effects and shock wave imprints are robust sig- 
nals for specific neutrino mixing scenarios that are un- 
likely to be mimicked by anything else. The information 
available from them is in the form of a combination of 



mass hierarchy and a range, so that complementary 
information from long baseline experiments is also re- 
quired. If 0i3 and mass hierarchy is already determined at 
terrestrial experiments, concrete information on the pri- 
mary neutrino fluxes will be obtained. On the other hand, 
if the burst is observed before the mixing parameters are 
measured, we shall have some advance idea about the 
neutrino mixing parameters expected, and that will guide 
our efforts towards the neutrino factory. 
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